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ABSTRACT 

In a previous work (Pichardo et al. 2005), we studied stable configurations for cir- 
cumstellar discs in eccentric binary systems. We searched for "invariant loops" : closed 
curves (analogous to stable periodic orbits in time- independent potentials) that change 
shape with the binary orbital phase, as test particles in them move under the influence 
of the binary potential. This approach allows us to identify stable configurations when 
pressure forces are unimportant, and dissipation acts only to prevent gas clouds from 
colliding with one another. We now extend this work to study the main geometrical 
properties of circumbinary discs. We have studied more than 100 cases with a range 
in eccentricity ^ e ^ 0.9, and mass ratio 0.1 ^ q ^ 0.9. Although gas dynamics may 
impose further restrictions, our study sets lower stable bounds for the size of the cen- 
tral hole in a simple and computationally cheap way, with a relation that depends on 
the eccentricity and mass ratio of the central binary. We extend our previous studies 
and focus on an important component of these systems: circumbinary discs. The radii 
for stable orbits that can host gas in circumbinary discs are sharply constrained as a 
function of the binary's eccentricity. The circumbinary disc configurations are almost 
circular, with eccentricity < 0.15, but if the mass ratio is unequal the disk is offset 
from the center of mass of the system. We compare our results with other models, and 
with observations of specific systems like GG Tauri A, UY Aurigae, HD 98800 B, and 
Fomalhaut, restricting the plausible parameters for the binary. 

Key words: circumstellar matter, discs - binary: stars. 



1 INTRODUCTION 

It is currently believed that fragmentation is the most prob- 
able mechanism for star formation, and the main product 
of fragmentation are multiple stellar systems with prefer- 
ence for wide eccentric binaries with separations ^ IQAU 
(Bonnell & Bastien 1992; Bate 1997; Bate & Bonnell 1997; 
Bodenheimer et al. 2000). Even in isolated stars, there is 
evidence that the majority of Sun-like stars formed in clus- 
ters (Carpenter 2000; Lada & Lada 2003), including the Sun 
(Looney, Tobin & Fields 2006). 

In the last decade the interest in binary systems has in- 
creased. This is in part because of the discovery that many 
T- Tauri and other pre-main sequence binary stars possess 
circumstellar and circumbinary discs as inferred from obser- 
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vations of excess radiation at infrared to millimeter wave- 
lengths, polarization, and both Balmer and forbidden emis- 
sion lines (Mathieu et al. 2000, Itoh et al. 2002, for a review 
see Mathieu 1994). On the other hand, recent observations 
of binary star systems, using the Spitzer Space Telescope, 
show evidence of debris discs in these environments (Trilling 
et al. 2007) and planets (Fischer et al. 2008). In their stud- 
ies they find that 60% of the observed close binary systems 
(separations smaller than 3 AU) have excess in their thermal 
emission, implying on-going collisions in their planetesimal 
regions. 

Over 150 extrasolar planets have been identified in sur- 
veys using the Doppler technique. Of the first 131 extra- 
solar planetary systems that have been confirmed, at least 
40 are in binary or multiple systems (for an up-to-date list 
see Haghighipour 2006). Approximately 30 of them are on 
S-type orbits (around one of the components: circumstel- 
lar discs) with wide stellar separations (between 250 and 
6500 AU), including at least 3 that orbit one member of a 
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triple star (Raghavan et al. 2006). Although most of these 
binaries are very wide, a few have separations smaller than 
20 AU (Els et al. 2001; Hatzes et al. 2003), challenging 
standard ideas of Jovian planet formation. Some interesting 
ideas try to explain the formation of Jovian planets within 
close binaries, but they could only explain few cases (Pfahl 
& Muterspaugh 2006), if more are discovered soon, these 
theories would not be sufficient. Although close binaries are 
not included in precise Doppler radial velocity search pro- 
grams because of their complex and varying spectra, at least 
one planet with a minimum of 2.5 Jupiter masses has been 
detected in a P-type orbit (around both components: cir- 
cumbinary discs), with a distance from the center of mass 
of 23 AU. The source is a radio pulsar binary comprised by 
a neutron star and a white dwarf in a 191 day stellar orbit 
(Lyne et al. 1988; Sigurdsson and Phinney, 1993; Sigurds- 
son et al. 2003). An example of accretion in P-type orbits 
about close binaries is given by Quintana & Lissauer (2006) 
who note the observation of the two small moons orbiting 
in nearly circular/planar orbits about the binary system 
Pluto-Charon (Weaver et al. 2006). Specifically regarding 
to circumbinary disc material, millimeter and mid-infrared 
excess emission has been detected around several spectro- 
scopic pre-main sequence binary star systems including GW 
Ori (Mathieu et al. 1995), UZ Tau E (Jensen et al. 1996), 
DQ Tau (Mathieu et al. 1997). 

In this work we have followed the same steps as in 
Pichardo, Sparke & Aguilar (2005, hereafter Paper I), where 
we opted for a simpler approach, analogous to using the 
structure of periodic orbits in a circular binary, to predict 
the gas flow. The path followed by a gas parcel in a stable 
disk around a star must not intersect itself, or the path of 
a neighboring parcel (unlike the case of planets, where the 
paths may cross) . In our work we follow Rudak & Paczynski 
(1981) and explore those non-crossing orbits of test particles 
that could be interpreted as gas particles in the low pressure 
regime, or as protoplanets or planets. An important issue in 
Celestial Mechanics is to determine the regions around a 
stellar binary system where accretion discs can form. Im- 
portant theoretical effort carried out to answer this ques- 
tion is reviewed in Paper I, where we studied circumstellar 
and circumbinary discs in binaries of arbitrary eccentricity 
and mass ratio. In this work we extend those studies, which 
were based on identifying families of stable invariant loops, 
a concept introduced by Maciejewski & Sparke (1997, 2000) 
in studies of nested galactic bars. We focus this time specif- 
ically on the geometry of circumbinary discs. We employ 
for this approach a test particle method probing the orbital 
structure of binaries of various eccentricities and mass ra- 
tios. 

In Section [5] we brieffy review the concept of an invari- 
ant loop, describe the method to solve the motion equations 
and the strategy used to find invariant loops. The geometry 
of the circumbinary discs including a fit for the inner radii 
of the circumbinary disc and a fit for the lopsidedness, are 
presented in Section |3] In Sections |4] and [5] we apply this 
study to compare with theoretical work and observations of 
some well known systems, respectively. Our conclusions are 
presented in section (6] 



2 THE METHOD AND NUMERICAL 
IMPLEMENTATION 

A more detailed description of the invariant loops method 
and its numerical implementation is given in Paper I (also in 
Maciejewski & Sparke -1997, 2000-). We give in this section 
a brief description. 

In the well studied circular 3-body problem, one known 
integral of motion is conserved: the Jacobi constant, defined 
in the rotating reference frame of the stars. Stable periodic 
orbits in this rotating system are defined and represent the 
"backbone" of the orbital structure. On the other hand, 
when the eccentricity is non-zero, there are no known inte- 
grals of motion to facilitate any analytical studies. However, 
the lack of a global integral of motion, does not preclude the 
existence of restrictions that apply to particular orbits. For 
motion in the plane of the binary, an additional integral of 
motion would confine an orbit to lie on a 1-dimensional curve 
every time the system comes back to the initial orbital phase. 
For example, if we look at the system every time the binary 
is at periastron, a particle following this orbit will land in a 
different spot but on the same 1-dimensional curve, which 
we call an invariant loop. In this manner invariant loops rep- 
resent the generalization of periodic orbits for periodically 
time-varying potentials. This means that an invariant loop 
is not a simple orbit but an ensemble of orbits that lie on a 
3-torus in this extended phase-space, but supported by an 
additional isolating integral of motion that forces the parti- 
cles to have a 1-D intersection with the orbital plane at a 
fixed binary phase. 

The equations of motion for the binary system are 
solved in terms of the eccentric anomaly ijj (Goldstein, Poole 
& Safko 2002, Section 3.7). We use units where the gravi- 
tational constant G, the binary semi-major axis a, and its 
total mass mi 4- m2 are set to unity so that, the binary 
period is 2-n, and its frequency oj = 1. The separation be- 
tween the stars at time t, measured from periastron where 
the azimuthal angle 9 = G, \s given by the radius r, 

r = a(l — ecostp) , (1) 
ujt = ~ e sin tp) , (2) 
cos 6^ = a(cos — e)/r. (3) 

The binary eccentricity, defined as e = a/1 — IP- jd? 
where a, and h are the semimajor and semiminor axes, and 
the mass ratio q = m-ilimx + m-2) are the only free param- 
eters. We use an Adams integrator (from the NAG fortran 
library) to follow the motion of a test particle moving in the 
orbital plane of the two stars. Kepler's equation ((2)| is solved 
with a tolerance of 10"''. 

The equations of motion of the test particle are solved 
in an inertial reference frame using Cartesian coordinates, 
with their origin at the centre of mass of the binary. All 
test particle trajectories are launched when the binary is at 
periastron, with the two components lying on the x-axis. 
The computation is halted if the particle runs away, moving 
further than 10 times the semimajor axis from the centre 
of mass, or if it comes within a distance of either star that 
results in a high number of force computations, in general 
due to close approaches to the stars. 

To find stable invariant loops, for which the phase space 
coordinates of our test particle, traces a 1-dimensional curve 
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on successive passes through periastron, we launch parti- 
cles from a chosen position along the a;-ajds joining the two 
stars at periastron, and examine the iterates in some two- 
dimensional subspace, such as the x — y plane. We plot the 
positions of the test particle at each complete binary pe- 
riod, and adjust the starting velocity Vy until the iterates 
converge on a one-dimensional curve. In practice, we look at 
the scatter along the radial direction for those iterates that 
lie within a sector that spans a small angle (5°) about the x 
axis when viewed from the centre of mass of the system. We 
adjust the launch velocity Vy until the radial scatter of pe- 
riastron positions of the test particle in the 5° sector drops 
within a threshold value. A value of 10~*a is used for cir- 
cumbinary loops and 10~^a for circumstellar loops. These 
values are consistent with the numerical errors in the orbit 
integration. For the majority of orbits, 10 points within the 
sector suffice and no more than 5 attempts are necessary to 
identify a given loop. While we are in a region of stable in- 
variant loops, the required launch velocity Vy is a continuous 
function of the starting point x. 

The numerical strategy we employ to solve the problem 
will allow us to find only stable invariant loops. Particles 
launched close to an unstable loop would diverge and the 
code would not be able to find this kind of loops. However, it 
is the stable orbits we are interested in. We do not calculate 
all the possible loops, but restrict our attention to those that 
are symmetric about the line joining the two stars when they 
are at periastron. When the binary orbit is circular, these 
are exactly the closed periodic orbits of a circumbinary disk. 
Although our figures show a set of discrete curves, invariant 
loops form continuous families in the same way as periodic 
orbits. We show only a few of the possible invariant loops, 
for clarity. 



3 GEOMETRICAL CHARACTERISTICS OF 
CIRCUMBINARY DISCS 

Discs in multiple stellar systems have attracted attention 
because of their high abundance, and also due to the in- 
teresting effects of the interaction on the disc morphology, 
better studied every day with the improvement of observa- 
tions. Of particular interest are the circumbinary discs that, 
because of their low density, are very difficult to observe. 
However, their importance arises from the possibility that 
these envelope discs might be feeding putative circumstellar 
discs, or harbor protoplanets, planets or any kind of debris. 

As in Section 3 of Paper I, the inner radius of the cir- 
cumbinary disk is set by the criterion that the stable loops 
exist and do not intersect each other. If the binary eccen- 
tricity is not small (e > 0.1) then the loops become unstable 
before they begin to intersect each other; we find no more 
one-dimensional curves, or the test particles fall toward the 
stars or go out of the system. 

To characterize the geometry of circumbinary discs, we 
use the coefficients of the Fourier expansion of the innermost 
loop: 



Ak = — ^ s(0i)cos(fc0i), 



Bk^ s[4>i)sin{k4,^) , (4) 

i = l 

where {si,(j)i) are the polar coordinates of A'' evenly spaced 
(in cf>) points along the innermost stable loop, measured from 
the binary center of mass. The modulus A\ + is used 
to determine the mean distance to the barycentre (fc — 0) 
and the lopsidedness {k — 1). 

More than one hundred simulations were included to 
calculate the fits we present in the next two subsections that 
provide the main geometrical characteristics of the discs. 
These simulations were performed in the eccentricity inter- 
val e — [0.0, 0.9], and mass ratio q = [0.1, 0.5] (equivalent to 
sample the whole range q — [0.1,0.9] because of the sym- 
metry in the definition of q). It is worth to mention that 
we have not included in the fits values for g < 0.1 since the 
behavior of the radius at those extreme values of q, changes 
abruptly and requires many more calculations. We will pro- 
duce fits for extreme q values (as is the case of planets) in 
a further paper. The loops technique, however, allows us to 
reach these extreme cases and we present in this paper an 
example of it (Section [5}. 

At any time, the invariant loops form closed curves, 
which deform as the binary follows its orbit, and return to 
their original shape when the binary returns to the same 
phase. But we found in Paper I that in practice even the cir- 
cumstellar loops do not deform strongly, and the circumbi- 
nary loops even less. Thus we measure the properties of the 
cirumbinary disks, and show their shapes in the figures, at 
the time of periastron passage. 



3.1 Lopsidedness of circumbinary discs 

An interesting effect produced by a high binary eccentricity 
combined with a large mass contrast, is the displacement of 
the geometric centre of the circumbinary disc (inner edge) 
with respect to the barycentre (See for example Figure [l|. 
This effect is a physical characteristic that could explain 
some observed asymmetries in discs (Duchene et al. 2004; 
Boden et al. 2005; Kalas et al. 2005). The displacement of 
the disc centre may affect calculations of the disc's inclina- 
tion (e.g. Itoh et al. 2002) since it is usual to link asymme- 
tries to inclination effects rather than to intrinsic asymme- 
tries in the geometrical centre of the discs. 

In theoretical work of dynamics in planetary systems, 
there have been several studies using linear analysis of per- 
turbations, which consider low mass ratios or low eccentrici- 
ties (e.g. Wyatt et al. 1999, Kuchner & Holman 2003; Deller 
& Maddison 2005). In these studies it is shown how the pres- 
ence of a second body (like a planet) in a system with an ec- 
centric orbit would impose a forced eccentricity on the orbits 
of the constituent dust particles, thus shifting the geometric 
centre of the disc away from the mass centre of the binary. In 
this manner, the dust in slightly eccentric orbits would glow 
more brightly when it approaches the pericentre. This could 
explain the asymmetry in the double-lobed feature in many 
systems (Holland et al. 1998; Koerner et al. 1998; Schneider 
et al. 1999; Telesco et al. 2000; Kalas et al. 2005; Freistetter 
2007). In the same direction, Dermott et al. (1999) find that 
if there is at least one massive perturber in the HR 4796 sys- 
tem that is on an eccentric orbit, then the system's secular 
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Figure 1. Circumbinary disc computed with invariant loops for 
a binary with q = 0.1 and e = 0.7, viewed at the moment of 
periastron marked with the two stars at {x, y) fs; (0, 0). The orbits 
of the stars are shown with the darker curves. The disc is not 
centred about the centre of mass of the system 



perturbations could cause the geometric centre of the disc 
to be offset from the star. 

Our technique of invariant loops allows the displace- 
ment of the disk to be calculated for arbitrary mass ratio 
and eccentricity. For mass ratios g 0.1 we have obtained a 
fit for the displacement of the circumbinary geometric centre 
with respect to the centre of mass of the binary (see Figure 



Rsh{e, q) = -Ciae' (0.5 - q) [q (1 - q)Y' 



(5) 



where a is the semimajor axis Ci = 3.7, 5 = 0.8, and r] = 
0.25, give our best fit to the calculated off-centre distance. 
The displacement in our calculations is always directed to 
the left of the centre of mass due to the position we chose to 
place the primary (left) and secondary (right) with respect 
to the centre of mass of the binary at their pericentre (see 
Figure [T| . 

We find the best fit to the proposed equation by search- 
ing the minimum residuals (some of them are given in Table 
[l]) produced by the square root of the sum of the squared 
differences between the proposed function, equation [5] and 
the computed displacements. The standard deviation com- 
paring the calculated data and the fit is 0.029a. The residual 
values in Table [T] are in units of the semimajor axis, a. 

The inner circumbinary rim is not exactly elliptical; 
close to resonances, the loops can even become slightly tri- 
angular. However, unless the eccentricity is very close to 
zero, the shape of the inner rim is close to circular and its 
maximum and minimum diameters lie almost perpendicu- 
lar. The eccentricity reaches a maximum value Cd ~ 0.15, 
almost independent of the binary eccentricity. 
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Figure 2. Lopsidedness measured by the quantity ij^j, (e, q)/e'', 
as a function of eccentricity e and mass ratio q, from equation [5] 
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Figure 3. Contour plot of the computed average inner radius 
RcbIo- for the circumbinary discs, from equation [6] 



3.2 Inner Radii (the "gap") 

In Paper I, we derived a simple relation for the size of cir- 
cumstellar discs as a function of binary mass ratio and ec- 
centricity. In this paper we extend the study to circumbinary 
discs and derive a relation now for the inner radius ( "gap" ) . 
Together, these radii define the region where no loops exist. 
We found that the change in radius of the circumbinary disc 
with the binary phase is very small; so here we calculate the 
inner radius of the circumbinary discs at one phase, when 
the stars are at their periastron. 

The best fit to our calculated inner radii is given by 



RcB'ye,q)-^C2a{\ + ae^) (g(l - g))" 



(6) 



where a is the semimajor axis, C2 = 1.93, a = 1.01, fi = 
0.32, and 7 = 0.043. The standard deviation obtained is 
0.09a. 

In Figure O we show a contour plot of the approxima- 
tion to the average radii for the circumbinary discs from 
equation |6] In Table [2] we show the difference between the 
relation |6] and the computed radii. 
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Table 1. Difference between the shift of the center of the circumbinary disk computed using invariant loops, and the fit (equation |5j , in 
units of the semimajor axis, a for some chosen pairs (e,q). 



e 




0.1 


0.2 


0.3 


0.4 


0.5 


0.00 




0.000 


0.010 


0.006 


0.010 


0.001 


0.20 




-0.005 


-0.025 


0.030 


0.009 


-0.004 


0.40 




0.012 


0.011 


-0.001 


0.009 


-0.001 


0.60 




-0.012 


-0.052 


0.012 


0.024 


0.000 


0.80 




0.001 


0.009 


0.081 


0.015 


0.000 



Table 2. Difference between the computed inner radii of the circumbinary disks and the fit (from relation [6ll , in units of the semimajor 
axis, a, for some chosen pairs (e,q). 



e q: 0.1 0.2 0.3 0.4 0.5 



0.00 


-0.051 


-0.065 


-0.005 


-0.034 


0.067 


0.20 


0.000 


-0.060 


-0.016 


0.107 


-0.204 


0.40 


0.009 


0.250 


0.203 


0.216 


-0.297 


0.60 


0.170 


-0.022 


0.106 


0.105 


-0.202 


0.80 


-0.040 


-0.120 


-0.003 


-0.062 


-0.392 



The ratio of equations [5] and [6] gives an estimate of the 
size of the displacement of the circumbinary disc geometric 
centre, as a fraction of the inner circumbinary (average) ra- 
dius. In Table Owe present the ratio between the calculated 
shifts and the corespondent average circumbinary radius. 



4 COMPARISON WITH THEORETICAL 
WORK 

In the problem of accretion and planet formation in binary 
systems, it is of paramount importance to determine the 
regions where a circumbinary disc can exist. Extensive liter- 
ature exists, but mostly devoted to the circular orbit binary 
case. Several studies have been done about the dependence 
on mass ratio and orbital eccentricity for the existence and 
characteristics of circumbinary discs. In particular, Quin- 
tana & Lissauer (2006) have simulated the late stages of ter- 
restrial planet formation within circumbinary discs in close 
binary systems with a wide range of orbital eccentricities. 
Quintana & Lissauer (2007) simulated the final stages of 
terrestrial planet formation in S- and P-type orbits within 
main-sequence binary star systems. 

In planetary studies a lot of work has been done (e.g. 
Holman & Wiegert 1999, hereafter HW99; Wyatt et al. 
1999). In the last reference, in particular, the authors ad- 
dress this issue by investigating the long-term stability of 
planetary orbits numerically using the eccentric restricted 
three-body problem. They launch circular prograde orbits 
in the vicinity of the stars and in the circumbinary region, 
looking for orbits that remain close to the stars for more 
than 10,000 binary periods. They provide a fit for the outer 
radii of the circumstellar discs, and for the inner part of the 
circumbinary disc, that depends only on the parameters of 
the binary (a, e, q). Their study has the advantage that it is 
able to examine a full range of eccentricity and mass ratios. 
The disadvantages are the lower precision of this technique, 
due to the fact that a disc can live much more than the fidu- 
cial 10,000 binary periods used by the authors to qualify an 



orbit as stable, and the fact that it is expensive computa- 
tionally. 

In our work we identify stable non-intersecting loops 
where gas may accumulate, a disc may develop and a planet 
may form. This condition plays the same role as searching 
for stable non-intersecting periodic orbits in the potential of 
a circular binary; in either case, the results could be modified 
by pressure or viscous forces. In this sense our search is more 
stringent than the study of HW99, who find orbits where a 
planet may survive for long times around a binary, but these 
orbits are permitted to intersect themselves or neighboring 
orbits. Still, a comparison of our results with those of HW99 
is relevant to gauge to what extent our different criteria 
result in similar constraints. 

We have calculated the difference between the fits made 
by HW99 for the circumstellar and circumbinary discs. In 
the case of circumstellar discs we find a good agreement with 
their results. For the circumbinary discs there are some dif- 
ferences. In the Figure|3]we show the results of HW99 (filled 
triangles) and the results obtained from invariant loops (con- 
tinuous lines) of the calculated average inner radii of cir- 
cumbinary discs, including the minimum and maximum dis- 
tance from the center of mass (open circles), vs the binary 
eccentricity, for two mass ratios, q = 0.1 (left), q = 0.3 
(right). We see in the figure that the fit by HW99 gives in 
general larger radii for the gap, especially at higher eccen- 
tricities of the central binary. For smaller eccentricities the 
radii provided for the fit of HW99 are almost the same or 
even smaller than the ones provided by the invariant loops. 
It is likely that as the binary becomes more eccentric, the 
phase space in which orbits can be trapped so that they must 
remain close to the stable circumbinary loops shrinks in vol- 
ume. That would make it less probable that the initially- 
circular orbits of HW99 would lie in that trapped region. 

5 APPLICATION TO OBSERVATIONS 

Lim & Takakuwa (2006) have already applied our results in 
Paper I to constrain the eccentricity of the binary in L1551 
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Table 3. Displacements calculated from invariant loops for some cases used for the fit in equation [S] given as a percentage of the 
corresponding circumbinary average radius. 



e 




0.1 


0.2 


0.3 


0.4 


0.5 


0.00 




0.000 


0.557 


0.320 


0.538 


0.051 


0.20 




8.291 


7.859 


3.785 


2.038 


0.149 


0.40 




12.340 


9.682 


7.177 


3.363 


0.035 


0.60 




15.217 


15.509 


9.235 


4.219 


0.000 


0.80 




20.179 


17.178 


9.577 


5.776 


0.000 




Figure 4. Comparison between the radii calculated by the fit of 
HW99 (filled triangles), and the fit with invariant loops that gives 
the mean radius (continuous lines) as eccentricity increases. We 
have also included the maximum and minimum radius from the 
centre of mass due to the shift of the disc (empty circles). Left 
panel: mass ratio q = 0.1, right: mass ratio q = 0.3 

IRS5 to e < 0.3, based on the sizes of the circumstellar and 
circumbinary disks. 

As a further application of our study, we have chosen 
four systems. The first represents the prototype of circumbi- 
nary discs: GG Tauri A. The second is UY Aurigae, the third 
is HD 98800 B, and the last is Fomalhaut. 

5.1 The Circumbinary Disc of GG Tauri A 
System 

GG Tan is a well known young multiple system. The sys- 
tem has two binary stars: GG Tau Aa/Ab, and GG Tan 
Ba/Bb. GG Tau A is an interesting binary since it possesses 
a circumbinary disk resolved in the millimeter wavelengths 
(Kawabe et al. 1993; Guilloteau et al. 1999). It has been ob- 
served at high resolution in the optical (Krist et al. 2002), 
and in near-infrared (Roddier et al. 1996; Silber et al. 2000; 
McCabe et al. 2002). 

The structure of the circumbinary disc seems to be char- 
acterized by an annulus with an inner radius between 180 
and 190 AU (Duchene et al. 2004; Guilloteau et al. 1999), 
and an outer radius extending up to 800 AU. The total 
mass of the circumbinary material {H2 + dust) is « 0.12 
Mq (Guilloteau et al. 1999). Itoh et al. (2002) derive an 
inclination of approximately 37°, assuming the orbit of the 
binary is coplanar to the circumbinary disc. 

The orbital characteristics of the central binary are still 
controversial. While Roddier et al. (1996), propose an ec- 
centric orbit in which the stars are located near periastron, 
at the same observation time, Krist et al. (2002) find they 



are close to apoastron in a highly eccentric orbit. McCabe 
et al. (2002) deduce an elliptical orbit e = 0.3 ± 0.2 with a 
semimajor axis of a = 35lg^ AU . The mass of each com- 
ponent of GG Tau A a and b, binary is obtained by White 
et al. (1999): 0.78±0.1 Mq, and 0.68±0.03 Mq respectively 
(mass ratio q ~ 0.47). 

Based on equation [6l we have calculated a band of pos- 
sible solutions for different semimajor axis and eccentrici- 
ties that could give an inner radius of 180 ± 18 AU that 
we present in Figure [5] In the same figure we locate the 
prediction by McCabe et al. (2002), whose error bar in the 
semimajor axis locates it inside our error zone (those with 
e = 0.3 and semimajor axis around 55 AU). In the same 
manner, some of the values derived by the predictions based 
in observations of Itoh et al. (2002), rest inside or close to 
this region, specifically those with semimajor axis a = 50 
AU , and eccentricities e = 0.4, 0.5. 

In figure[6]we present two possible configuration for GG 
Tau constructed with invariant loops. The first plot is us- 
ing the values of eccentricity and semimajor axis reported 
by McCabe et al. (2002) (triangle on Figure [5} . Notice that 
this configuration results in a significantly reduced inner gap 
which is not compatible with the result reported by Duchene 
et al. (2004) and Guilloteau et al. (1999). In the second 
panel, we present a configuration that is compatible with 
these two references {Rcb ~ 180 AU) and still within the 
error box of McCabe around the preferred values of (e = 0.5, 
a = 53 AU). In both cases we have assumed a mass ratio of 
q = 0.47. Although this parameter has little importance in 
determining the circumbinary disc inner radius, it is impor- 
tant in fixing its lopsidedness. 

The circumbinary disc seems to be shifted away from 
the center of mass of the binary (Duchene et al. 2004) by 
~ 0".16 or 22 AU for a distance of 140 pc. We predict that 
the shift of the disc geometrical centre with respect to the 
barycentre should be small since the stars have almost the 
same mass. As shown in Figure O we find that the shift pro- 
duced by the binary system on the circumbinary disc, would 
be less than 4 AU. We calculate that a mass ratio of g < 0.3 
would be necessary to produce both the observed shift and 
the observed inner radius of the circumbinary disk. This is 
unlikely because the effective temperature and luminosity of 
the two stars are very similar (White at al.l999). 

Because the observational parameters determined by 
McCabe et al. (2002) and Itoh et al. (2002) imply that the 
circumbinary disc should extend much closer to the binary 
and the lopsidedness should be much smaller than is pre- 
dicted by considering orbits around the binary GG Tau A, 
alternative theories have arisen. An interesting idea is the 
possibility that the smaller GG Tau B binary (10" away 
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Figure 5. Possible configurations for GG Tau A with our model 
using the observed inner radius of the circumbinary disc and an 
error from observations of ±18 AU (shadowed region), assuming 
the disc is only sculpted by the binary. We have also included 
some observational and model predictions for the configuration. 
The filled squares are the predictions in the work of Itoh et al. 
2002. The filled triangle is the prediction of McCabe et al. 2002 
including the error bars. Finally we show five possible curves (con- 
tinuous lines labeled with different Rgh) with the calculated shifts 
(equation [Sjl , for the given q = 0.47, and the corresponding pair 
(a, e). 



and with masses 0.12 ± 0.02 and 0.044 ± 0.006 Mq for GG 
Tau B a and b respectively) has something to do with the 
sculpting of the circumbinary disc (Beust & Deutrey 2006). 
Unfortunately these authors find that, although there are 
possible arrangements to explain the outer edge of the disc, 
the binary B is not massive enough, and not close enough 
to explain the wide inner radius in this way. 




-WO 



a 

x(AS..I) 



HID BOO 

q-0.A7. a -0.5, 0-63 



5.2 The Binary System UY Aur 

UY Aur is a binary system of classical T Tauri stars 
(Duchene et al. 1999). It is located in the Taurus- Aurigae 
star forming region at an approximate distance of 140 pc 
(Elias 1978). The projected separation on the sky between 
UY Aur primary and secondary (A and B), is 120 AU (Close 
et al. 1998). The spectral types of UY Aur A and B are 
estimated to be MO and M2.5, respectively (Hartigan & 
Kenyon 2003). Assuming a circular orbit, the binary pe- 
riod is ~ 1640 ± 90 yr and the total mass of the binary is 
~ 1.73 ±0.29 Mq (Hioki et al. 2007). 

A circumbinary disc around the UY Aur binary was 
detected by near-infrared, polarimetric and millimeter inter- 
ferometric ^^CO emission observations (Duvert et al. 1998). 
Its inner radius is about 520 AU (Hioki et al. 2007). These 
authors deproject the circumbinary disk assuming that its 
inner edge is circular and the binary is coplanar with it. 



Figure 6. a) Circumbinary disc computed with invariant loops 
for the central value of eccentricity (e = 0.32) and scmimajor 
axis (a = 35 AU) derived by McCabe et al. (2002) for a mass 
ratio q = 0.47. b) Same as a) but for a plausible invariant loops 
solution for eccentricity (0.5) and semimajor axis (a = 53 AU) 
that reproduces the inner edge observed, of 180 AU, assuming 
that the only factor sculpting the inner edge of the circumbinary 
disc is the main binary, GG Tau A 



to get a binary separation of 167 AU and an inclination of 
42° ± 3°. The disc seems to be not uniform showing clumpy 
structure, circumstellar material inside the inner cavity and 
an arm-like structure, probably created by accretion from 
the outer region of the disc or stellar encounters (Hioki et al. 
2007). 

We calculated stable invariant loops in a system with 
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Figure 7. Possible configurations for UY Aurigae with our fit 
for the circumbinary radius for a mass ratio q = 0.25, using the 
observed circumbinary disc inner radius of 520 AU and an error 
of ±20 AU (shadowed region), assuming the disc is only sculpted 
by the binary. We show nine possible curves (continuous lines 
labeled with different Rsh) with the calculated shifts (equation 
[Sjl, for the given q = 0.25, and the corresponding pair (a, e). The 
darker curves are the lines a(l + e) = 167 AU and a(l — e) = 167 
AU 



the orbital parameters proposed for the binary by Hioki et al. 
(2007), a binary separation of 167 AU and e = 0. We have 
set q = 0.5, but Figure |3] shows that the edge of the CB disk 
does not change significantly for any q > 0.1. This model 
would give an inner radius for the circumbinary disc of no 
more than 340 AU , far less than the observed 520 AU for 
this system. 

We have proceeded by constructing a family of solu- 
tions as in the case of GG Tau A (Section 15. l|l . This is, we 
provide a family of solutions which reproduce the approx- 
imate inner radius of the circumbinary disc. We construct 
two different possible configurations, with different mass ra- 
tios, q = [0.25, 0.4] (Figures[7]and|8l), which bracket q = 0.36 
as given by Hartigan & Kenyon (2003). Although the inner 
radius of the circumbinary disc is almost independent of 
the mass ratio (for q ^ 0.1), the shift of the circumbinary 
disc increases as the mass ratio becomes more unequal. The 
present separation is 167 AU , which must lie between the 
periastron and apastron separations. This restricts e, a to 
points above the curves given by a(l — e) and a(l + e) in 
Figures [7] and |8l Thus the binary cannot be circular, but 
must have e > 0.1 and for a mass ratio of g = 0.36 (Harti- 
gan & Kenyon 2003), the center of the disk should be offset 
from the mass center of the stars in the direction towards the 
center of the secondary star's orbit by Rsh ^ 0.05a. Since 
the apparent separation has remained constant since 1944 
(Hioki et al. 2007), the system is unlikely to be very close 
to periastron, implying a larger eccentricity. 
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Figure 8. Same as Figure [7] but for q = 0.4 

5.3 The Binary System HD 98800 B 

HD 98800 (HIP 55505, TWA 4A) is a hierarchical quadru- 
ple star system in the TW Hya association. The separation 
between A and B components is approximately 0.8" on a 
north-south line (Prato et al. 2001). Torres et al. (1995) 
find that both visual components are themselves spectro- 
scopic binaries. Tokovinin (1999) derived orbital parameters 
a = 62 AU, e = 0.5 for the orbit of HD 98800 B around HD 
98800 A. 

The system HD 98800 B has excess flux in the mid- 
infrared, which was interpreted by Soderblom et al. (1998) 
and by Prato et al. (2001) as a circumbinary dust disk. Their 
estimates of the inner radius of the disk were 1.5 to 2 AU. 
New observations of the stellar binary were reported by Bo- 
den et al. (2005). They estimated visual and physical orbits 
of the HD 98800 B subsystem with interferometric observa- 
tions combined with astrometric measurements by the Hub- 
ble Space Telescope Fine Guidance Sensors. The orbital and 
physical parameters obtained in that work are given in Table 

HI 

Using the parameters inferred by Boden et al. (2005), 
Akeson et al. (2007) used the results in our Paper I to argue 
that the average inner radius of the circumbinary disk should 
be at RcB ~ 3.4 AU. Thus the central hole in the cirumbi- 
nary disk should be significantly larger than suggested by 
Soderblom et al. (1998) and Prato et al. (2001). Treating 
both binaries (A and B) as point masses, each with the 
combined mass of its two stars, q « 0.5 and from Table 
1 in Paper I we predict that the maximum outer radius of 
the circumbinary disc of B is about a tenth of the semi- 
major axis of the orbit of A and B, or 6.2 AU. Akeson et al. 
compare these results with a dynamical simulation of the 
system, in which the two stars of binary B are followed sep- 
arately while binary A is modeled as a single star, and the 
circumbinary disk of B is represented by test particles. After 
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Table 4. Approximate Orbital Parameters of HD 98800 B from Boden et al. (2005). 



q 0.45 

e 0.78 

a 0.98 AU 

Period (days) 314 
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Figure 9. Stable invariant loops making up the circumbinary 
disc of HD98800 B, using the parameters of Boden et al. (2005), 
given in Table|3] The orbits of the stars are shown with the darker 
curves 



the equivalent of 1 Myr, the inner edge of the disk was at 
about 3 AU and the outer edge at 10 AU , in approximate 
agreement with predictions from our invariant loops. 

We have constructed the corresponding circumbinary 
disc for HD98800 B (Figure [gj. Using equations [S] and [5] we 
find that the disc center is shifted with respect to the centre 
of mass of the system by Rsh ~ —0.1 AU . The shift is small 
because the two stars have nearly equal mass. 



5.4 The Fomalhaut Dust Belt 

Fomalhaut (HD 216956 or a Pisces Austrinus) is a bright 
nearby A3 V star of 2 Mq, an age of 200±100 Myr (Bar- 
rado & Navascues 1998) at a distance of 7.7 pc. It shows a 
dust ("debris") ring around it between 133 and 158 AU from 
the central star (Aumann 1985; Guillet 1985; Dent et al. 
2000; Holland et al. 1998, 2003) with an approximate mass 
between 50 and 100 Earth mass (Kalas et al. 1995). This 
structure represents one of the best observed extrasolar ana- 
log to our Kuiper belt. Fomalhaut's ring has an inclination 
of 24° away from edge on. The disc presents an asymmetry 
in the brightness with the southern side nearer the star than 
the opposite side (Stapelfeldt et al. 2004; Marsh et al. 2005). 
The deprojected asymmetry (off-centred) of « 15 AU , the 
sharp inner cut at 133 AU , and the slight eccentricity of the 



disc Bd ~ 0.1 has been studied recently by Kalas, Graham & 
Clampin 2005 and Quillen (2006) who have proposed that 
all the characteristics of this system can be explained by the 
presence of a planet just interior to the ring inner edge. 

Several theories to form the sharp inner edge of the 
disc in Fomalhaut have been proposed by different studies, 
the most accepted one until this moment is the presence of 
inner planets. The solution for Fomalhaut is degenerate in 
the sense that several combinations of the main parameters 
(mass ratios between the central star and the planet, eccen- 
tricities and semimajor axes) could reproduce the observed 
values for the inner radius of the dust disk Rcb, and the 
displacement Rsh of its center from the star's position. We 
explore here three possible set of parameters given by other 
authors. 

In Table [5] we show the selected parameters for which 
we have calculated the invariant loops that would make up 
a circumbinary disc, and their corresponding references. 

Under the assumption that a secondary body (a giant 
planet) is the only source of asymmetry of the ring in Fo- 
malhaut we have calculated sets of invariant loops corre- 
sponding to both models in Table [5] In the upper panel of 
Figure [TO] are given the results with invariant loops for the 
approximation of Deller & Maddison (2005). From Figure [TOl 
and Table [6] we show that the inferred value of eccentricity 
from Deller & Maddison (2005) is considerably larger than 
it should be to obtain the observed inner radius of Fomal- 
haut's disk, and the shift, if the mass of the planet is 2 Mj. 
In the lower panel of Figure [TD] and second row of Table |S1 
we show a solution for a planet mass of 2 Mj for which Rcb 
and Rsh are closer to the observed values for Fomalhaut. We 
have also computed the eccentricity Cd of the inner rim of 
the circumbinary disk presented in the last column of Table 

El 

Quillen (2006) proposed solutions with a less-massive 
planet, using the formulae of Wyatt et al. (1999) to calcu- 
late the expected eccentricity. For two extreme values of the 
range given in that work, this is, the mass of Saturn and 
the mass of Neptune, given in Table [S] we have calculated 
the invariant loops to find the average and the shift radii. 
The values given by Quillen (2006) are close to the results 
that we derive from the invariant loops, as can be seen in 
Table [6] This shows that in the limit of small eccentricity, 
calculations based on invariant loops agree with those from 
earlier methods which are valid only in those regimes. 



6 CONCLUSIONS 

We have extended the studies started in Paper I to a more 
detailed analysis of circumbinary discs in eccentric binary 
systems from the geometrical and physical point of view. 
The disks are defined by a family of stable invariant loops, 
the analogs to stable periodic orbits around a circular bi- 
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Table 5. Orbital parameters proposed for Fomalhaut. The first column is the mass, the second is the planetary mass (in terms of the 
mass of Jupiter, Saturn or Neptune), eccentricity e and semimajor axis a are in the third and fourth columns, the employed technique 
and reference is given in the last two columns 





Mpi 


e 


a 


Technique 


Reference 


(Mq) 






(AU) 






2.3 


2 Mj 


0.4 


59 


N-body 


Deller & Maddison (2005) 


2 


1 Ms to 1 M;v 


0.1 


119 


Secular perturbations 


Quillen (2006) 



Table 6. Disc characteristics from invariant loops for Fomalhaut assuming the proposed masses of DcUer & Maddison (2005) and Quillen 
(2006) for the planet, and our best approximation to the solution taking the proposed masses by these authors. The first column shows 
the assumed mass for Fomalhaut, the second is the proposed mass for the planet around Fomalhaut (J=Jupiter, S=Saturn, N=Neptune), 
the third and fourth are the eccentricity and the semimajor axis of the planetary orbit, the fifth and sixth are the shift radii calculated 
with invariant loops, the seventh column is the value of the circumbinary inner rim eccentricity (ej), and the last column indicates the 
author of the paper we have taken the parameters from or our best approximation to the parameters that approach the most to the 
observed values of Fomalhaut. 





Mpi 


e 


a 


Rsh 


RCB 


fid 


Work 


(Mo) 






{AU) 


{AU) 


{AU) 






2.3 


2 Mj 


0.4 


59 


36 


155 




Deller & Maddison (2005) 


2.3 


2 Mj 


0.16 


61.5 


16 


133 


0.15 


Our best approx. 


2 


Mjv 


0.1 


119 


12.5 


142.5 




Quillen(2006) 


2 


Mjv 


0.13 


90.5 


14.5 


136 


0.13 


Our best approx. 


2 


Ms 


0.1 


119 


12.5 


148 




Quillen(2006) 


2 


Ms 


0.13 


106.4 


14 


135 


0.13 


Our best approx. 



nary, which do not cross each other or themselves. Thus 
they define paths that can be followed by clouds of gas, 
which dissipate energy when they run into each other. Just 
as with a binary in circular orbit, the circumbinary disk is 
truncated at its inner edge when there are no longer any sta- 
ble non-crossing orbits for the gas to follow. We have used 
this property of the invariant loops to define the inner edge 
of the circumbinary disk. 

We showed already in Paper I that the inner radius of a 
circumbinary disk depends strongly on eccentricity, opening 
wider gaps for higher eccentricities. The size of the inner 
hole depends only slightly on the mass ratio. The geometric 
centre of the circumbinary disc is off-centre with respect to 
the centre of mass of the binary system. The disc is closer to 
being symmetrical around the whole orbit of the secondary 
star. 

Here, we have explored the range of parameters to quan- 
tify both the off-centering of the circumbinary disc with re- 
spect to the centre of mass of the system, and the average 
inner radius of circumbinary discs, as a function of mass ra- 
tio and eccentricity. We compare our results with the work of 
HW99 who searched for initially-circular orbits that survive 
more than 10,000 periods of the binary. When the eccentric- 
ity is small this procedure gives similar results to ours, but 
at larger eccentricity HW99 find fewer stable orbits close to 
the binary, and hence larger inner gaps. This could be re- 
lated to the fact that the larger the eccentricity, the smaller 
the available phase space of orbits that are trapped so that 
they must remain close to a stable invariant loop. 

If the properties of a circumbinary disc arc observable, it 
is possible to constrain the binary system properties by com- 
paring the predictions based on invariant loops with what is 
observed. Likewise, if the orbital parameters of a binary are 



known, the geometry of the circumbinary and circumstellar 
regions permitted for stable orbits are readily obtained. 

For the well known circumbinary disc of the binary sys- 
tem GG Tau A, we use the inner radius of the circumbinary 
disk to restrict the possibilities for the binary parameters. 

Since two stars have nearly equal masses, we would expect 
the ring to be nearly symmetrical about the center of mass of 
the stars. The observed offset is substantial and much larger 
than can be explained by orbital dynamics; effects such as 
the finite ring thickness may be important (e.g. Duchene 
et al. 2004). 

In the case of UY Aurigae, we show that the observed 
radius of the circumbinary disk requires that the binary orbit 
be noncircular, with eccentricity e > 0.1. The center of the 
disk should be offset from the mass center of the stars in the 
direction towards the center of the secondary star's orbit by 
Rsh > 0.05o. 

We have modeled the system HD 98800 B with the pa- 
rameters given by Boden et al. (2005): the shift of the disc 
with respect to the centre of mass of the system should be 
about 0.1 AU. 

Although the fits we provide here are valid for values 
of g ^ 0.1, our technique is also applicable to extreme cases 
of g ^ 0.001. For the disk around Fomalhaut, we have used 
invariant loops to propose plausible solutions for the orbital 
parameters of a planet that explains the morphology of the 
debris disc of this system. 

We reach similar results to Quillen (2006), but for the 
larger planetary mass proposed by Deller & Maddison (2005) 
we find that the planet's orbit must be less eccentric. 
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Figure 10. Circumbinary disc of Fomaliiaut. Upper panel: in- 
variant loops disc using the parameters (e = 0.4, a = 59 AU) in 
Table [5] by Deller & Maddison (2005). Lower panel: our solution 
for the same planetary mass but with (e = 0.16, a = 61.5 AU). 
The trajectory of the planet is shown with the darker curves. 
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